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Abstract. The European Pulsar Timing Array (EPTA) network is a col- 
laboration between the five largest radio telescopes in Europe aiming to study 
the astrophysics of millisecond pulsars and to detect cosmological gravitational 
waves in the nano-Hertz regime. The advantages and techniques of handling the 
multi-telescope datasets of a number of sources will be presented. In addition, 
the results of the EPTA timing analysis of the pulsar-white dwarf binary PSR 
J1012+5307 will be reported. Specifically, the measurements for the first time 
for this system, of the parallax, the variation of the projected semi-major axis 
and of the orbital period. Finally, the derived stringent, theory independent 
limits on alternative theories of gravity, with the use of this ideal laboratory for 
strong- field gravity tests, will be presented. 



1. Introduction 

Millisecond pulsars (MSPs) are extremely stable and accurate cosmic clocks due 
to their short spin periods and the high degree of rotational stability. The tech- 
nique of pulsar timing is the regular monitoring of the rotation of these objects 
by tracking the times of arrival (TOAs) of the radio pulses. Frequent measure- 
ments of pulse arrival times of pulsars provide a powerful tool to determine their 
spin and astrometric parameters to a very high degree of precision. 

High precision timing of radio millisecond pulsars, most of which are mem- 
bers of binary systems, can be the key to a number of unanswered problems in 
fundamental physics and astronomy, ranging from stellar evolution to tests of 
gravitational physics in the strong field regime and the detection of a cosmolog- 
ical gravitational wave background. 



2. The European Pulsar Timing Array 



The European Pulsar Timing Array (EPTA) network (Table [I]) is a collabo- 
ration between the five institutes (ASTRON, JBO, INAF, MPIfR and Nancay 
observatory) operating the largest radio telescopes in Europe. It is consisting 
of the Effelsberg 100m radiotelescope of the Max-Planck-Institute for Radioas- 
tronomy (MPIfR) in Germany, the 76m Lovell radiotelescope of the Manchester 
University, at Jodrell Bank, UK, the 94m equivalent Westerbork Synthesis Radio 
Telescope (WSRT) of ASTRON in the Netherlands, the 94m equivalent Nangay 
decimetric radio telescope (NRT) in France and soon the 64m Sardinia Radio 
Telescope (SRT) in Italy. 
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The EPTA is using the available telescopes for high-precision timing in a 
coordinated way, in schedules and source lists. This results, in larger and denser 
datasets in multiple frequencies. In addition, exchange of data, people and 
knowledge between the working groups is ordinarily happening, resolving swiftly 
any systematic telescope problems or other issues. The main aim of the EPTA 
is to increase the precision and quality of pulsar timing measurements, to study 
the astrophysics of millisecond pulsars and to detect cosmological gravitational 
waves from coalescent massive binary black holes in the nano-Hertz regime. 



Table 1. The European Pulsar Timing Array collaboration. On the first 
column the institutes and telescopes and on the second the people participat- 

ing in the EPTA. 

Institute-Telescope People 



MPIfR-Effelsberg 

U. Manchester /Jodrell Bank 
-Lovell 

Nancay observatory-NRT 
INAF-SRT 

ASTRON/U. Amsterdam/ 
U. Leiden-Westerbork 



M. Kramer, A. Jessner, N. Wex, 

P. Freire, D. Champion, K. Lazaridis 

B. Stappers, A. Lyne, C. Jordan, 

G. Janssen, M. Purver, S. Sanidas 

I. Cognard, G. Theureau, 

R. Ferdman, G. Desvignes 

A. Possenti, M. Burgay, M. Pilia 

J. Hessels, Y. Levin, 

R. van Haasteren 



Currently the EPTA is collaborating with the Parkes Pulsar Timing Array 
(PPTA) in Australia and the nanoGrav in the USA, forming the International 
Pulsar Timing Array (IPTA), a global pulsar timing array with full sky coverage 
wishing to lead the way in the detection and study of gravitational waves. 

2.1. Effelsberg millisecond pulsars 

In the current work, several millisecond pulsars have been chosen from the Effels- 
berg sou rce list. On the archival data of those, new calibration techniques were 
applied (jLazaridid 120091 ) . i.e. ephemeris update, noise removal from individual 



channels, new time corrections, summation of individual scans and creation of 



new synthetic templates ( Kramer et al.|[l994l ). This re-calibration resulted in an 



improved timing accuracy making the pulsars ideal candidates in the effort of 
directly detecting gravitational waves in the nano-Hertz regime. In total fifteen 
sources, observed monthly with the Effelsberg radio telescope at at least two dif- 
ferent frequencies, have been chosen. In Table [5] the selected sources are shown 
with the current post-fit rms of the model fit achieved, in comparison with the 
post-fit rms without the improved calibration procedures. In all the cases the 
improvement is vast, varying from two times to two orders of magnitude. In the 
last column the post-fit rms with the use of only the 1.4GHz data is shown. 

At least twelve of the analyzed ms pulsars can be good candidates for the 
EPTA efforts in detecting the stochastic background of gravitational radiation. 
Work is in progress for combination of the datasets from all the EPTA telescopes 
for all of these sources, which will improve the current measurements and take 
us closer to the desired detection limit. Until reaching this point the combined 
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Table 2. The fifteen analyzed millisecond pulsars with the post-fit rms with 
and without calibration. The TOAs are until September 2008 (MJD 54720) 
and all the available Effelsberg frequencies have been used. On the last column 
only the 1.4 GHz TOAs have been used. 



source 


Post-fit rms 
before (/xsj 


Fost-nt rms 
after (//s) 


Post-fit rms 

-t A /~\ TT / \ 

1.4 LxHz {fj,s) 


PSR J0030+0451 


a 


5.7 


3.8 


PSR J0218+4232 


51.5 


9.6 


9.1 


PSR J0613-0200 


28.1 


2.7 


2.6 


PSR J0621+1002 


23.8 


6.8 


6.5 


PSR J0751+1807 


15.5 


4.9 


4.3 


PSR J1012+5307 


3.5 


2.7 


2.6 


PSR J1022+1001 


16.2 


3.7 


3.1 


PSR J1024-0719 


~1500 


13.4 


2.7 


PSR J1518+4904 


23.5 


18.9 


19.0 


PSR J1623-2631 


~2100 


3.9 


3.9 


PSR J1640+2224 


15.2 


1.7 


1.5 


PSR J1643-1224 


34.4 


3.9 


3.8 


PSR J1744-1134 


1.7 


0.6 


0.6 


PSR J2051-0827 


49.9 


48.7 


48.8 


PSR J2145-0750 


4.1 


2.6 


2.5 



a : The source was regularly observed from 2008. 



datasets are being used for doing high precision timing analysis of the individual 
systems, as it is shown in the following section for PSR J1012+5307. 



3. PSR J1012+5307 



3.1. Introduction 

PSR J1012+5307 is a 5.3 m s pulsar in a low ec centricity binary system with 
orbital period of Pb =14. 5 h (IN icastr o et alJll995h and a low mass he lium white 
dwarf (WD) companion ( Lorimer et al.lll995al ). ICallanan et al.l ( 19981 ) compared 
the measured optical luminosity of the WD to the value expected from WD 
models and calculated a distance of d = 840±90 pc. In addition they measured, a 
radial velocity component of 44±8kms _1 relative to the solar system barycenter 
(SSB), and the mass ratio of the pulsar and its companion q = m p /m c = 10.5 ± 
0.5. Finally they derived a companion mass of m c = 0.16 ± 0.02 M Q , a pulsar 
mas s of m v = 1.64 ± 0. 22 Mq and an orbital inclination angle of i = 52° ± 4' 



Lange et al.l ( 200ll ) presented a complete precision timing analysis of PSR 



J1012+5307 using 4 years of timing data from Effelsberg and 7 years from Lovell 
radio telescope. They derived the spin, astrometric and binary parameters for 
the system and they discussed the prospects of future measurements of Post- 
Keplerian parameters (PK) which can contribute to the derivation of stringent 
limits on alternative gravity theories. 
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In this work PSR J1012+5307 has been revisited with seven more years of 
high-precision timing data and combined datasets from the EPTA telescopes. 
The data have been analyzed using the timing software TEMPCQ. 

3.2. Results 

Combining almost 3000 TOAs, at five different frequencies from the four tele- 
scopes now in use by the EPTA, we have been able to improve on the timing 
solution and on all the astrometric, spin and binary parameters of this system, 
see Figure [TJ 




MJD 



Figure 1. Post-fit timing residuals using datasets from Effelsberg, Lovell, 
Nancay and Westerbork. 

For the first time a parallax ir = 1.2 ± 0.3 mas has been measured for 
PSR J1012+5307. This corresponds to a distance of d = 822 ± 178 pc which is 
consistent with the d = 840 ± 90 pc measured from the optical observations. In 
addition, the proper motion measurements have been improved by an order of 
magnitude, with a total proper motion of \i% = 25.735 ± 0.019 mas yr" 1 yielding 
a total transverse velocity of vt = 102.0 ± 9.8 kms -1 . Using the radial velocity 
from the optical measurements, the space velocity of the system v space = 111.4± 
9.5 kms" 1 has been calculated, being consistent and almost three times more 
precise than the previous value. 

An upper limit of e < 8.4 x 10" 7 (95 per cent C.L.) has been obtained, being 
consistent with the evolutionary scenario of spin-up of PSR J1012+5307 through 
mass transfer from the companion, while it is in the red giant phase. A change 
in the projected semi-major axis of x a b s = 2.3(8) x 10~ 15 has been measured 
in the current analysis, which has been used to constrain for the first time the 
position angle of the ascending node. 

As predicted bv lLange et al. ( 200ll ) a significant measurement of the change 
in the orbital period of the system, P& = 5.0(1.4) x 10 -14 , has been obtained for 
the first time. This is caused by the Doppler correction, which is the combined 



1 http: / /www.atnf.csiro.au /research / pulsar /tempo/ 
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effect of the proper motion of the system ( Shklovskii 19701 ) and a correction 
term for the Galactic acceleration and by a contribution due to the quadrupole 
term of the gravitational wave emission, as predicted by general relativity (GR). 
After subtracting these two contributions from our measured value, the excess 
value of P£ xc = (—0.4 ±1.6) x 10 -14 confirms the validity of GR for one more 
millisecond pulsar binary system. 

All the terms mentioned above are the ones expected to contribute by us- 
ing GR as our theory of gravity. However, there are alternative theories of 
gravity, that violate the strong equivalence principle (SEP) and predict extra 
contributions to the observed orbital period variation. One is the dipole term 
of the gravitational wave emission, which results from the difference in gravita- 
tional binding energy of the two bodies of a binary system. Thus, the case of 
PSR J1012+5307, where there is a pulsar- WD system, is ideal for testing the 
strength of such emission. One finds for small-eccentricity pulsar- WD systems, 
where the sensitivity s (related to the gravitational self-energy of a body) of 

the WD is much smaller than the one of the pulsar, Pt, %v ° e = — Air 2 -^r^ kds p 2 , 

( Willi l200ll ) where T & = 4.9255 jj,s and \i is the reduced mass. s p is the sensi- 
tivity of the pulsar and kd refers to the dipole self-gravitational contribution, 
which takes different values for different theories of gravity (zero for GR). 

Another term is predicted by the variation of the locally measured grav- 
i tational constant a s the universe ex pands, P^ = —2 ^ [1 — (1 + ^) s p ] Pf,, 
(jDamour et alJll988l : lNordtvedtlll990h where M is the total mass of the system. 
It has been shown that there is no need to add these extra contributions to 
explain the variations of the orbital period, however the excess value has been 
used to set limits for a wide class of alternative theories of gravity. 

PSR J1012+5307 is an ideal lab for constraining the dipole radiation term 
because the WD nature of the companion is affirmed optically, the mass es- 
timates are free of any explicit strong-field effects and the mass of the pul- 
sar is rather high, which is important in the case of strong field effects that 
occur only above a certain criti cal mass, like the spontaneous scalarisation 
(|Damour fe Esposito-Eares3l993h . Thus, by usi ng the G/G = (4±9)x 10" 13 yr" 1 
limit from the Lunar Laser Ranging (LLR) ( Williams et al.l r2004f ) the G/G 
contribution has been calculated and subtracted from our excess value in or- 
der to finally obtain an improved generic limit for the dipole contribution of 
k d = (0.2 ± 2.4) x 10~ 3 (95 per cent C.L.). 

A generic test for G cannot be done with a single binary pulsar, since general 
theories that predict a variation of the gravitational constant typically also pre- 
dict the existence of dipole radiation. This degeneracy has been broken here in a 
joint analysis of PSR J1012+5307 and PSR J0437-4715 (jVerbiest et al.ll2008l ). 
two binary pulsar- WD systems with tight limits for Pf, and different orbital pe- 
riods. By applying equation ^ = -2% [1 - (1 + ff|) s p ] - 4vr 2 ^ K D s 2 p 

to both binary pulsars, and solving in a Monte-Carlo simulation (Figure [2]) this 
set of two equations, stringent and generic limits based purely on pulsar data 
and in the strong field regime have been obtained. With a 95 per cent C.L., 

% = (-0.7 ±3.3) x 10~ 12 yr" 1 and k d = (0.3 ±2.5) x 10~ 3 . In the future, more 
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accurate measurements of Pf, and distance of the two pulsars and WDs could 
constrain even more our derived limits. 
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Figure 2. Contour plots of the one and two a confidence regions on G/G 
and kd jointly. The elongation of the regions reflects the correlation due to 
the mutual dependence of the two pulsar- WD systems in this combined test. 

A more de t ailed work on this study of PSR J 1012+5307 can be found in 
Lazaridis et al.1 (120091 ) . 
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